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ABSTRACT 



Context. With certain assumptions, radiative energy transport can be modelled by the diffusion approximation. In this case, the 
Rosseland mean opacity coefficient characterises the interaction between radiation and matter. The opacity data are usually available 
in pre-tabulated form, and in the generation of the data one assumes a distinct heavy element mixture, which is usually a scaled solar 
(-H ' one. Therefore, presently available data is unable to cover the full parameter range of some astrophysical problems, in which the 

chemical composition of the medium being considered varies. 
I ■ Aims. We attempt to produce low temperature opacity data incorporating the effects of varied abundances of the elements carbon and 

O ' nitrogen. For our temperature range of interest, molecules represent the dominant opacity source. Our dataset covers a wide metallicity 

range and is meant to provide important input data for stellar evolution models and other applications. 
C/3 , Methods. We conduct chemical equilibrium calculations to evaluate the partial pressures of neutral atoms, ions, and molecules. 

Based on a large dataset containing atomic line and continuum data and, most importantly, a plethora of molecular lines, we calculate 
Rosseland mean opacity coefficients not only for a number of different metallicities, but also for varied abundances of the isotopes '^C 
and '"'N at each metallicity. The molecular data comprise the main opacity sources for either an oxygen-rich or carbon-rich chemistry. 
We tabulate the opacity coefficients as a function of temperature and, basically, density. 

Results. Due to the special role of the CO molecule, within a certain chemistry regime an alteration to the carbon abundance causes 
considerable changes in the Rosseland opacity. The transition from a scaled solar (i. e. oxygen-rich) mixture to a carbon-rich regime 
results in opacities that can, at low temperatures, differ by orders of magnitude from to the initial situation. The reason is that the mean 
opacity in either case is due to different molecular absorbers. Variations in the abundance of nitrogen have less pronounced effects 
but, nevertheless, cannot be neglected. 

Conclusions. In typical astrophysical applications, it is indispensable to take into account opacity variations due to chemistry changes. 
In this respect, the new data is superior to previous compilations, but is, however, still subject to uncertainties. 
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, 1. Introduction thermal conductivity, a basic en ergy transpor t proce ss in the in- 

5_( ■ ert C-O core, is available from ICassisi et al.l ( l2007l) . Tabulated 

. Radiative energy transport within a medium plays a role in a temperature opacity data (ranging from a few thousand to a 
wealth of astrophysical scenarios. When the photon inean free j^^^^^^^^^ ^^^-^^ ^^^^-^^ ^^^^^^^ ^ ^j^^ Qp^L coUab- 
path IS short compared with the typical scale height of the oration dlglesias& BLogerd 119961) and the Opacity Project (OP, 
medium, irrespective of the photon wavelength, this problem giii3a|2005D. Both groups provide their data publicly through 
can be modelled by the diffusion approximation. By assuming interface^ where in general single element abundances 
local thermodynamic equilibrium (LTE), the definition of a har- ^^^-^^ Qp^L also allows to produce so-called Type II 
monic meari opacity coefficient - the so-called Rosseland mean ^^^j^^. ^^^^ ^ ^^^^^^^ ^^^^^^^ ^j^^ abundances 
d Rosselandl 1 19251) - emerges naturally from the evaluation of ^^^^^^^^ ^ enhanced by adding con- 
ffie integrated radiative flux. Rosseland stated that the deriva- ^^^^^ ^^^^ ^^^^^-^^^ ^^^^ ^^^^ ^ ^^.^^ hov^cvcr, no 
tion of the analytical expression for this characteristic function j^olecular contributions to the Rosseland mean are taken into 
of the medium regarding its power of transmitting radiation en ^^^^^^^ ^^^^ Qp^^ and OP databases. Molecules be- 
bloc IS not difficult, yet the actual calculation is. This is be- ^^^^ dominant opacity source at temperatures lower than 
cause one usuaUy has to account in detail for a wide variety of ^^^^^ 4000-5000 K, which occur in the outer envelope and at- 
absorbers. Typical applications of Rosseland mean opacities ai-e j^osphere of a stai- on one of the giant branches - either the 
stellai- structure and evolution calculations, which span orders of ^-^^^ ^^^^^^ ^^^^^ _^ protoplanetary ac- 
magnitude in temperature and density and require the inclusion ^^^^-^^ ^-^^^ ^^^^ j^^^^ temperatures (below approximately 
of numerous absorption and scattering processes. As an exam- ^^qq ^^^^ responsible for the bulk of opacity, 
pie, we consider a low mass star in the late phase of its evolu- 
tion on the Asymptotic Giant Branch (AGB). Data for electron 

' OPAL: http : //www-phys . llnl . gov/Research/OPAL/ 
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One of the first papers giving a detailed discussion about 
low temperature mean opacities (and also summarising ear- 
lier efforts on this topic) was Alexander et al. (1983). This 
wo rk evolved further and result ed in the extensive database 
of Alexander & Fergusonl (1 19941) and the updated version of 
iFerguson et alj ( 20051) - henceforth AF94 and F05, respectively 
- that has become a standard for low temperature opacities in the 
past few years. These data are based on scaled solar metal com- 
positions (plus some tables for enhanced alpha element abun- 
dances). Returning to our AGB star example, the problem is 
that the chemical composition in the envelope of such an object 
varies. Products of the ongoing nucleosynthesis in the stellar in- 
terior are brought to the surface by a series of mixing events. The 
entire m echanism is dubbed Third Dredge-up (TDU, for a re- 
view see lBusso et al ][T999.) . The main burning products dredged 
up are freshly synthesised carbon and elements produced by 
the slow neutron capture process. Alterations to the element 
mixture result in significantly different opacity coefficients. The 
corresponding implica tions for the ste l lar evo lution calculations 
were emphasised by Alexander et al.l (Il983b . The authors pro- 
vide some examples of how the Rosseland mean opacity changes 
when the number ratio of carbon to oxygen atoms (C/0) varies. 
The important role of the C/O ratio is due to the fact that of all 
the abundant molec ules, carbon m onoxide (CO) has the high- 
est binding energy (lRussel]|[l934<) . The partial pressure of CO 
hardly varies with the C/O ratio in a plasma at constant temper- 
ature and pressure. In a chemical equilibrium situation, the less 
abundant species of carbon and oxygen are almost completely 
locked in the CO molecule. The remaining fraction of the more 
numerous atoms is free to form other molecules. The absolute 
number of free atoms (i. e. not bound in CO) of either oxygen or 
carbon primarily determines the magnitude and characteristics 
of the Rosseland mean opacity since CO contributes on average 
less than other molecular species (such as H2O, TiO, C2H2, and 
HCN). 

Although the above facts have been in principle known for a 
long time, they have not been accounted for in stellar evolution 
models, due to the lack of adequate opacity data. From the tabu- 
lated low temperature opacities that can be found in the literature 
(other than AF94 and F05, e. g. iN euforge 1993; Semen ov et al 
2003HWeiss et alj2006l:[Sharp & Bun'ows.2007;.Freedman et al 



20081 partly monochromatic and focused on s pecial objects 
such a s brown dwarfs and extrasolar planets), only lHelling et alJ 
(120001) considered a case in w hich C/O - 1.8 i n simulating 
winds of carbon-rich AGB stars. Potter et al. I (l200l investigated 
single element abundance variations and the influence on low 
and high temperature opacities exemplarily. Subsequently, the 
sensitivity of stellar evolution models to these changes was anal- 
ysed and both effective temperatures and stellar lifetimes depend 
clearly on opacity changes due to abundance variations. 

A dramatic improvement in evaluating molecular opacities 
with a varying amount of carbon was achieved by Marigo (2002) 
who illustrated how a correct treatment of opacity can affect stel- 
lar evolution tracks by using analytical AGB star models. The 
principle outcomes were that observations of carbon stars can 
be reproduced more accurately than before and the models ap- 
pear to imply that the carbon star phase is shortened considerably 
with a consequent reduction in the stellar yields. The opacities 
were estimated by chemical equilibrium calculations and ana- 
lytical fits for molecular contributions to the Rosseland mean, 
which were combined by a simple linear summation. This was 
the main problem of this approach. The non-linear nature of the 
Rosseland mean by definition renders the method of adding up 
opacity contributions a fragile approximation. Moreover, the de- 



rived molecular contributions are not unique because they cannot 
be determined in a non-ambiguous way. The Rosseland mean 
emphasises transparent spectral regions, and the gaps between 
absorption lines from a specific molecule can be filled by other 
species depending on the chemistry. Thus, this interplay cru- 
cially influences the total opacity. Finally, not all relevant molec- 
ular absorbers have been taken into account in the work de- 
scribed above. 

From this overview, we conclude that the stellar evolution 
community require a complete, homogeneous database contain- 
ing Rosseland mean opacity coefficients with varied abundances 
of carbon. The aim of this work is to fill this gap. We also in- 
cluded abundance variations in nitrogen, since in low and in- 
termediate mass stars mechanisms might activate the CN cycle 
and trigger the reconversion of carbon isotopes into nitrogen. 
The corresponding mechanism acting in stars of mass lower than 
4M0 is dubbed Cool Bottom Process (Wasser burg et al.|[T995l) 
which is a slow, deep circulation at the bottom of the convec- 
tive envelope. Its origin is still unknown, although a number 
of driving mechanisms have been proposed (e.g. rotationally 
induced instabilities, magnetically induced c irculation, gravity 
waves or thermohaline mixing - we refer to iBusso et al.l 120071 
for an overview). In more massive AGB stars (M > AMp), th e 
analogous process is known as Hot Bottom Burning (Ilbenll973h . 

A preliminary version of the opacity data presented in the 
following was applied sti c cessfu ll y by X^ ristallo et al. (200^ 
2008) and iLebzelter et all (l2008l) . ICristalio et al.1 (12007) con- 
sidered a low metallicity (Z = 1 x 10"*) stellar model with 
M = 2M0. By using the new opacity coefficients, the model 
was able to reproduce observational data more accurately than 
earlier models in terms of physical properties (such as for in- 
stance the effective temperature) and the abundance pattern of 
the heavy elements. Information th at is complementary to th e 
data considered here can be found in lLederer & Aringeil (l2008l) . 

The paper is structured as follows. In Sect.|2l we describe the 
tools used to generate our opacity tables. Section [3] summarises 
all data adopted in our calculations, i. e. abundances as well as 
atomic and molecular opacity data. We describe and justify the 
design of our tables in Sect. |4] We discuss the results in detail in 
Sect. |5] before concluding and providing a perspective on future 
work in Sect. |6l 



2. Tool and method 

To generate the da ta presented in this work, we used the COMA 
code developed bv lAringe3(l2000l) . For a description of improve- 
ments since the initial version, we refer to a forthcoming paper of 
Aringer et al. (in preparation). Assuming local thermodynamic 
equilibrium, the program solves for the ionis ation and chemical 
equilibrium (using the method of lTsuiilll973h at a given temper- 
ature and density (or pressure) combination for a set of atomic 
abundances. From the resulting partial pressures for the neutral 
atoms, ions, and molecules, the continuous and line opacity is 
calculated at the desired wavelengths using the data listed in 
Sect. [3] The main purpose of the COMA code was originally 
to provide monochromatic absorption coefficients for dynamical 
model atmospheres of cool giants (Hofner et al. 2003). However, 
it has been used in a wide range of applications, for example in 
the process of calculating low and high r esolution spectra (e. g. 
lAringer et alj|200a ILebzelter et al.ll20p 8l respectively) and line 
profile variations "("Nowotnv et ani2005h . 

Once the monochromatic opacities Ky(T,p) = K(v,T,p) are 
known {T stands for the temperature and p for the density), the 
calculation of a mean opacity coefficient such as the Rosseland 
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mean is straightforward. One has to perform a weighted integra- 
tion of Ky over the relevant frequency range. For the Rosseland 
mean /cr - kr{T,p), the relation is given by 



Table 1. Continuous opacity sources 



Ion and process 



Reference 



]__ n 1 dByjT) ^^ j r dBy{J)_^^ 
kr Jo Ky dT I Jo dT 



(1) 



where the weighting function is the partial derivative of the 
Planck function with respect to the temperature. The main sub- 
ject of this paper is to study not only the dependence of kr on the 
thermodynamic quantities T and p, but in addition the chemical 
composition. 

In practice, the integration over v in Eq.[T]must be performed 
at a predefined discrete frequency (or waveleng th) grid. We use 
a grid that is based on the one described by |j0rgensen et al.l 
( Il992h . but we extend the wavelength range to have boundaries 
at 200,000 cm-i (500 A = 0.05 yum) and 200 cm"' (50//m). This 
results in a total number of 5645 wavelength point s at which we 
calculate the opacity. iHelling & J0rgensenl (Il998b proposed to 
use a number of opacity sampling points for the accurate mod- 
elling of atmospheres of late-type stars that is roughly four times 
larger than the value used here. The same is true for the number 
of wavelength points in F05, who use more than 24,000 points. 
However, the error introduced when using a lower resolution is 
generally small compared with other uncertainty sources (see 
Sect. 15. 3t , and a smaller number of grid points has the advan- 
tage of a reduced computing time. We note that the opacity sam- 
pling technique remains - regardless of the precise number of 
points - a statistical method. To arrive at a realistic and com- 
plete description of the spectral energy distribution, one would 
need a far higher resolution (R ^ 200,000). In any case, the grid 
is sufliciently dense for a rectangle rule to be sufficient in car- 
rying out the wavelength integration. This can be justified by 
comparing the numerically obtained value of the normalisation 
factor on the right-hand side of Eq. [1] with its analytical value. 
The formal integration limits in the definition of the Rosseland 
mean have to be replaced by cut-off wavelengths. These are de- 
termined by the weighting function dBy{T)/dT that constrains 
the relevant wavelength range. Like the Planck function itself, 
the maximum of its derivative shifts to higher wavelengths with 
decreasing temperature and vice versa. At the upper wavelength 
limit adopted here (50;um) and for the lowest temperature at 
which we generate data (about 1600 K), the weighting function 
has decreased by more than 1/100 relative to its maximum value. 
Accordingly, at the high temperature edge the weighting func- 
tion at the same wavelength has dropped to almost 1/10,000 of 
its maximum value. Since we do not include grain opacity in our 
calculations that would require going to even higher wavelengths 
(as in F05), we definitely cover the relevant spectral range for 
the calculation of kr within the adopted parameter range. The 
decline in the weighting function towards lower wavelengths (or 
higher frequencies) is far steeper, so that the above argument is 
also fulfilled at the low-wavelength cut-off. 

3. Data sources 

In the following, we briefly summarise the sources of the data 
entering the opacity calculations. The basic ingredient in this 
procedure is the relative amount of elements contained in the 
mixture for which we would like to know the opacity. In this 
work, we chose to use the set of recommended values for solar 
element abundances compiled by Lodders (2003), which imply 
a solar C/O ratio of 0.501. The values f or the elements C, N , 
and O are close to the values given by iGrevesse et al.l (I2007D : 



Hib-f andf-f 
H- b-f 
H- f-f 

H+H (Quasihydrogen) 
H^ f-f 
H7 f-f 

Ci, Mgi, Ah, Si I, He I f-f 
He- f-f 

continuous scattering 
Rayleigh scattering from H i 
Rayleigh scattering from H2 



Karzas & Latter (1961) 
Doughty et al. (1966) 
Doughty & Fr aser(1966) 
Doyle (1968) 
Mihalas (1965) 
Someryille(1964) 
Peach (1970) 
Somerville (1965) 
Mihalas (1978) 

Dalgarno, quoted by Gingerich (1964) 
Dalgamo & WilUams (1962) 



using their abundances results in (C/O)© = 0.537. The authors 
derive these values from a 3D hydrodynamic model of the solar 
atmosphere, a technique that caused a downward revision of the 
solar CNO abundances in re cent years. However, these values 
are still disputed. For instance. ICenteno & Socas-Navarrol(l2008h 
used spectro-polarimetric methods to argue for an oxygen abun- 
dance that is higher than claimed by Grevesse et al. (20 03) and 
closer to previously accepted values (e. g. Greves se & Sauvall 
1998. with (C/0)o = ( ).490) that agree with values derived from 
helioseismology (e. g. lBasu & Antiall2008l) . 

However, the data presented depend to first order only on 
the relative amount of carbon and oxygen, which does not differ 
significantly from that for the various abundance sets mentioned 
above. Therefore and due to the fact that C/O is a variable quan- 
tity, the current tables can serve as an approxim ation for appli- 
cations that use abundances other than lLoddersI (I2003D . until we 
generate further data. 



3.1. Continuous opacity 

The routines for the calculation of the continuum opacity in 
COMA are based on an earlier version of the MARCS code 
(lJ0rgenserietal.lll992h. Th e latest MARCS release was de- 
scribed by Gustafsson et al.l ('2008). We adopt the format of their 
Table 1 to ease comparison with their updated set of continu- 
ous opacity sources. The data that we use (and Ust in Table [1]) 
is not as extensive. However, the most relevant sources are in- 
cluded, and in the low temperature region of the presented tables 
in particular the molecular contribution to kr dominates over the 
continuum by several orders of magnitude. 



3.2. Atomic lines 

Atomic line data are taken from the VALD database 
jKupka et al ]l2000lfL where we use updated version 2 data (from 
January 2008) here. For the atomic lines, we adopt full Voigt pro- 
files derived from the damping constants listed in VALD. The 
only exception is hydrogen for which we use an interpolation in 
tabulated line profiles from! Stehl3 (Il995h. The atomic partition 
functions are taken from the work of'irwin (198 1), although the 
data for boron has been modified (Gorfer 2005). The number of 
atomic lines included is 16,059,201. Split into their respective 
ionisation stages the numbers are: I 4,028,995, II 5,347 ,990, 
and III 6,682,216. For the ionisation energies, we refer to IStiftI 
(2000) since we use the same reference data quoted there in the 
current work. 



VALD: http : // ams . astro . univie . ac . at/ " vald/ 
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Table 2. Molecular line data 



Molecule 


Thermodynamic 


Number of lines 


Line 




data" 




data * 


CO 


1 


131859 


7 


CH 


2 


229 134 


8 


C2 


1 


360 887 


9 


SiO 


2 


85 788 


10 


CN 


1 


2 533 040 


8 


TiO 


1 


22 724 670 


11 


H2O 


3 


27 693 367 


12 


HCN/HNC 


4 


33 454021 


13 


OH 


1 


36 601 


14 


VO 


1 


3 171552 


15 


CO2 


2 


1032 266 


16 


SO2 


5 


29 559 


17 


HF 


1 


462 


18 


HCl 


1 


447 


17 


ZrO 


1 


16391 195 


19 


YO 


1 


975 


20 


FeH 


6 


116 300 


6 


CrH 


1 


13 824 


21 


C1YI2 




opacity sampling 


22 






opacity sampling 


22 



° (1) ISauval & Tatuiiil 11984 (2) iRossi etalJ Il985l (3) 
Vidle r & Tennyson' J2000l ~(4') iBarber et alj l200l TsT llrwirJ 1 19881. 

(6)|Dulick et al. 200l 

(7) Goorvitch & Chackerian 1994, (8) J0rsensen 1997, (9) 
Ouerci et al . 1974, (10) Lanahoff & Bauschlicher 1993, (11) Schwenke 
19981 (12)lBarberetal. 2006, (13) IHarri s et al. 2006, (14) Schwen^ 
1997i (15)lAlvarez&Plez. ■ 19981. (16 ) .Rothman et al. 1991 (17) 
Rothman et alj|2005l. (18)lTii3pindl2007l (19)|Plezll2007r(20)lLittletoj 



Kotnman et at JI2UU3I, ( 1 S ) 1 1 ii)pmai2UU /l il^) H^tezlUUU II ( ZU ) 
2007L (21) lBauschlicher et alJl200ll. (22) |j0rgensen et alJll98S 



3.3. Molecular data 

In the calculation of low temperature opacities, molecules play 
a critical role. We list the data that we used to calculate the ta- 
bles containing Rosseland mean opacity coefficients in Table |2] 
References to the thermodynamic data (i. e. the partition func- 
tion) and the line data for each molecule (displayed in the first 
column) are as indexed in the columns two and four, respec- 
tively. The number of lines entering the calculation is given 
in the third column (the original lists contained more lines). 
For some molecules, there is more than one line list available. 
The line lists that we use were selected during other projects 
(cited at the beginning of Sect. |2]l proven to deliver viable re- 
sults. In the case of the OH line list, the measured da ta from the 
HITR AN database contains more lines than that from lSchwenkd 
( 11997 ). We performed some test calculations after completing 
the database with OH data from HITRAN. The change in our 
results was, however, hardly perceptible. As long as the overall 
opacity distribution is reproduced reasonably, the line positions 
do not have to be precisely correct in calculating kr (as opposed 
to applications relying on high-resolution spectral synthesis). 

For the molecules C2 and CN, we intr oduced some modi- 
fications to the line data. |j0rgensenl (Il997h suggested a scaling 
of the gf v alues for C2 in a certain wavelength region (details 
are given in lLoidl et alj|200il) . We followed this suggestion but 
also investigated the effect on a-r of not applying this scaling 
(see Sect. 15.3b . The modifications to the CN line list are not 
crucial to the calculation of the Rosseland mean. We corrected 
the line positions of approximately 18,000 lines using measure- 
ments listed in the catalogue of Davis et al. (2005). The data 
for the molecules C3 and C2H2 from J0rgenseni (il997D are only 



available in the form of an opacity sampling. For the genera- 
tion of this data, a microturbulent velocity of 2.5 km s"' was 
adopted. As for the calculation of the molecular line opacity, 
we assume Doppler profiles in each case, since there is little in- 
formation about the damping constants of molecular transitions. 
Additionally, for the species dominating the overall opacity and 
in regions close to bandheads, even the wings of the strongest 
lines will contribute less to the opacity than the Doppler cores of 
the numerous overlapping neighbour lines. 

The chemical equilibrium constants use d in calculati ng of 
the molecule partial pressures are t hose f rom lTsuiil (Il973h with 
updates described in 'Helling et al. (1996*). In the case of C3, we 
utilised data published by Irwin (1981 ). All molecules (a total 
number of 314 species) enter into the equation of state. 



4. Table design 

We tabulate the logarithm of the Rosseland mean opacity 
log /(-R [cm^ g"'] as a function of log T[K], the logarithm of the 
gas temperature, and log/? [gcm^^K^^ 10'^], where R = p/T^. 
p and represents the density in units of g cm"^ and the tem- 
perature in millions of Kelvin, respectively. The ranges cov- 
ered are 3.2 < log T [K] < 4.05 with a step size of 0.05, and 
-7.0 < \ogR < 1.0 with a step size of 0.5. The low tempera- 
ture cut-off was set to be the temperature at which grains may 
become the major opacity source (cf. e. g. F05). However, typ- 
ical AGB stellar evolution models do not attain such low tem- 
peratures. We would like to emp hasise that du st formation is 
usually not an equilibrium process dGail & Sedlmavr.198 8), and 
thus an a priori tabulation of grain opacities can be affected by 
large uncertainties. At the highest value of log T, the contribu- 
tion of molecules to the mean opacity vanishes and a smooth 
transition to high temperature opacity data is possible (see Sect. 
15.11 ). Data are available for 14 different values of Z (cf. Table[3]) 
that depicts the total mass fraction of all elements heavier than 
helium. The metallicity spans a range fro m Z - \ x 10"^ to the 
Magellanic clouds metallicity range ( van den Berg h 2000a bJ), 
where the grid is a bit denser in terms of Z, and the solar metallic- 
ity to super-solar metallicity (Z = 0.04). We calculated tables for 
three different mass fractions of hydrogen {X € {0.5,0.7,0.8}). 
Thus, the data cover the supposedly prior field of application, 
namely the outer layers of an AGB star, where such values of X 
occur Initially, we calculated a master table for each metallic- 
ity, in which the metal composition was linearly scaled from the 
abundances given by Lodders (2003) to arrive at the respective Z 
value. From these abundances, we then enhanced the mass frac- 
tions of '^C and '"^N in steps that depended on the initial metal- 
licity (see Table[3]). Since this produced an increase in the overall 
metallicity Z, we followed the OPAL approach and reduced the 
mass fraction of '^He to fulfil X+Y + Z - 1 . 

The number of enhancement factors is constrained by a 
trade-off between numerical accuracy and computational costs 
in the generation and application of the tables. Due to the special 
role of carbon in conjunction with oxygen (see Introduction), the 
number of enhancement factors is higher than for nitrogenQ We 
calculated opacities for 7 different mass fractions of '^C. The 
starting point was the mass fraction that results from scaling all 
element abundances to the metallicity under consideration. All 
other carbon mass fractions resulted from multiplying this mass 
fraction by factors chosen as follows. For each metallicity, we 



H ere we deviate from the original design used in ICristallo et al.l 
j2007b and lLederer& AringeJ j2CI08l) . where we adopted 5 enhance- 
ment factors for both carbon and nitrogen . 
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Metallicity (Z) '^C enhancement factors '''N enhancement factors 



0.04 


1.0 1.2 


1.5 


1.8 


2.0 


2.2 


3.0 






1.0 


1.2 


1.5 


0.02 


1.0 


1.5 


1.8 


2.0 


2.2 


3.0 


5.0 




1.0 


1.5 


2.0 


0.01 


1.0 


1.5 


1.8 


2.0 


2.2 


5.0 


10.0 




1.0 


1.5 


3.0 


0.008 


1.0 




1.8 


2.0 


2.2 


5.0 


10.0 


20.0 


1.0 


2.0 


4.0 


0.006 


1.0 




1.8 


2.0 


2.2 


5.0 


10.0 


30.0 


1.0 


2.5 


5.0 


0.005 


1.0 




1.8 


2.0 


2.2 


5.0 


10.0 


30.0 


1.0 


2.5 


5.0 


0.004 


1.0 




1.8 


2.0 


2.2 


5.0 


10.0 


30.0 


1.0 


2.5 


5.0 


0.003 


1.0 




1.8 


2.0 


2.2 


6.0 


20.0 


70.0 


1.0 


3.0 
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use the factors 1.8, 2.0, and 2.2. The C/O ratio emerging from 
the adopted scaled solar abundances was about 0.5. Thus, multi- 
plying the initial X('^C) by 2.0, we inferred that C/O - 1, where 
the molecular opacity in general reaches a minimum at low tem- 
peratures. As one can see for example in Fig. [8] the molecular ab- 
sorption increased sharply on both sides of this minimum, while, 
towards much higher and lower C/O ratios, there was some type 
of saturation. To resolve this sharp turnaround, we included the 
factors 1.8 and 2.2 (corresponding to C/O - 0.9 and 1.1, respec- 
tively). In stellar evolution models of low metallicity AGB stars 
this is probably of minor importance because the transition to a 
carbon star is usually rapid and occurs within a few dredge-up 
episodes (e. g. iMarigo & Girardiir2007h . The respective highest 
enhancement factor is related to the initial metallicity. From the 
work of Abia et al. (2003) and referenc es therein and the first ap- 
pHcation of our data bv iCristallo et al.l fcOO?! [20081) . we derived 
information about the final carbon abundances reached in AGB 
stars. These determined the maximum enhancement factors in 
our tables for the metallicities under consideration (for instance 
Z - Ix 10"', 1 X 10"\ 1 X 10"**). At all other metalHcities, we de- 
rived the highest factors using a roughly linear relation between 
log Z and log /cmax (where we denote the enhancement factor as 
/). The remaining two factors for carbon are distributed almost 
equispaced on a logarithmic scale between f - 2 and /c.max 
at the lower metallicities. For high Z with low final enhance- 
ments, we shifted factors from the carbon-rich regime to the re- 
gion in which C/O < 1. For '^N, we introduced two additional 
factors beyond the initial abundance. The expected final over- 
abundance of nitrogen is much lower than for carbon, and we 
set /c,max//N.max = 10 for the lowcst valucs of Z, decreasing this 
value for increasing Z. The intermediate enhancement factor for 
nitrogen was set to bisect the logarithmic interval approximately 
between the two other factors. An overview of all enhancement 
factors at each metallicity is given in Table[3] 

We have not yet considered varied alpha element abun- 
dances, although one expects an incre ased value of [O/Fe] at 
low metallicities (e. g. iHill et alj|2000l and references therein). 
The reason is the sharp increase in the number of tables, if one 
retains the previously outlined data structure. In its current ver- 
sion, the database contains 3 x 7 x 3 = 63 tables per metallicity. 
Varying the abundance of oxygen also influences the C/O ratio, 
which is the decisive quantity for the molecular opacity at low 
temperatures. This in turn requires alterations to the enhance- 
ment factors of carbon, since one wishes to retain at least one 
point where C/O = 1, even for an enhanced oxygen abundance. 
Establishing a scheme with a minimal number of enhancement 



factors for three elements (or element groups) is therefore not 
straightforward if one is attempting to retain as much informa- 
tion as possible with respect to the role of the C/O ratio. In place 
of enhancement factors, one could add constant amounts of one 
element (group) in terms of a mass fraction, as completed in the 
OPAL Type II tables. Only the application of the data in its cur- 
rent form will provide us with information about the feasibility 
of our approach and whether the data should be arranged in a 
different way. Future work will be dedicated to these questions. 

5. Results and discussion 

The calculation of Rosseland mean opacities for scaled solar 
metal mixtures has been discussed extensively in many papers 
(see Introduction for citations). Since we find good agreement 
with data from other groups (see Sect. 15. It , we only briefly re- 
state the main points of this procedure. In Fig. [T] we provide 
an overview of the contents of the database using examples. We 
refer to individual panels in the following paragraphs. Generally 
speaking, molecules cause the mean opacity to vary dramatically 
as a function of temperature in a similar way for each hydrogen 
mass fraction X and metallicity Z (left panels of Fig. [T]i that we 
include in our database. Beyond a temperature corresponding to 
log T - 3.6 to 3.7 (depending on Z and log/?), the contribution 
of molecules to kr vanishes and only continuous sources and 
atomic lines block the radiation field. 

AF94 provided an in-depth discussion about which type of 
opacity is dominant in different regions of the parameter space. 
They presented a detailed treatment of the main molecular opac- 
ity sources, which were in this case water (H2O) - accounting for 
the large bump at low temperatures - and titanium oxide (TiO). 
They also provided information about the monochromatic ab- 
sorption coefficients of these molecules. Beside H2O and TiO, 
there is a number of other molecules that deliver non-negligible 
cont ributions to the Ro sseland opacity in the oxygen-rich case 
(cf. iLederer & Aringerli200 8). At higher temperatures, CN and 
CO contribute to kr, and VO, OH, and SiO (ordered by decreas- 
ing importance) should also be taken into account. Calculations 
based on these 7 molecules result in opacity coefficients that 
are accurate to about 10 per cent compared with the full dataset 
when the metal mixture is oxygen-rich. The further inclusion of 
CrH and YO reduces this error to below 3 per cent on average. 
The temperature ranges in which different molecules contribute 
to kr indicated in the bottom right panel of Fig. [1] are estimated 
by omitting the respective molecules in the calculation of kr and 
by checking where the change in this quantity exceeds 5 per cent 
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Fig. 1. Contents of the opacity database in some showcases. Top left panel: Rosseland mean opacities at constant log 7? for Z = 0.02 
as a function of log T for different values of X. The qualitative behaviour is fairly independent of X and Z and varies smoothly with 
log/? at fixed abundances of carbon and nitrogen (the latter is not enhanced in any of the panels shown in this figure). Bottom left 
panel: The molecular opacity decreases when lowering the metallicity, but the structure with a bump at low temperatures due to 
the molecular contribution to the opacity is conserved. Top right panel: For the case where carbon is enhanced the molecules also 
produce high opacities at low log T, although the shape of the curve differs noticeably from the standard case. Bottom right panel: 
The reason for the different structure in the opacities is that different molecules contribute to the opacity in the carbon-rich case 
(white boxes), in contrast to the oxygen-rich case (dark-grey boxes). The CO molecule contributes in both cases as well as CN, 
although at different orders of magnitude in each case. The extension of the boxes provide an indication of the temperature where 
the sources contribute but does not contain information about the order of magnitude of the contribution. The total contribution of 
molecules and atoms is assessed by leaving out these opacity sources in the calculations. All curves have been smoothed using cubic 
splines. See text for details. 



with respect to the complete dataset. However, the limits derived 
from this criterion vary with the value of log R under considera- 
tion. We use log/? = -1.5, which is typical of the envelopes of 
AGB stars (Sergio Cristallo, private communication). 

An increase in the carbon content of the metal mixture while 
the amount of oxygen remains at its original value (i. e. an in- 
crease in the C/O ratio) has distinct effects on the Rosseland 
mean opacity. We refer to the Introduction and the previous sec- 
tion for a description of the respective mechanism. In the transi- 
tion from the oxygen-rich regime (Fig.lT] top left) to the carbon- 
rich regime (Fig.IT] top right), one can distinguish between two 
cases. At lower temperatures (below log T - 3.4; we refer in the 
following to the case shown in Fig. |2l i.e. log/? - -3.0), the 
opacity first decreases, due to the above described property of 
CO, which causes the following mechanism. As C/O increases 
from its initial value of about 0.5 and approaches 1, more oxygen 



becomes bound in CO and fewer oxygen atoms are free to form 
molecules with a large overall absorption such as H2O. Close 
to C/O = 1 (not necessarily at an exact equal amount of C and 
O), the opacity reaches a minimum as the partial pressures of 
oxygen-bearing molecules drop substantially, while the abun- 
dances of carbon-bearing molecules only begin to rise to sig- 
nificant levels. As the amount of carbon continues to increase, 
thus incrementing C/O beyond 1, the opacity increases due to 
the formation of polyatomic carbon-bearing molecules such as 
C2H2 or HCN. These polyatomic molecules are obviously most 
relevant at lower temperatures. In addition, C3 and C2 produce a 
bump in the opacity at high carbon enrichments. 

The situation at higher temperatures (up to log T - 3.1) is, 
however, different. Besides CO, the contribution of which re- 
mains almost constant, only CN and C2 are relevant opacity 
sources (cf. Fig. [1] bottom right panel), while other molecules 
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Fig. 2. Rosseland opacity changing with the carbon content in the metal mixture as a function of temperature. We show curves 
for two metallicities, Z - 0.02 and the lowest metallicity in this database Z - 0.00001, at a value of log/? = -3.0. The full line 
represents the solar scaled metal mixture with C/O - 0.5. An increase in the carbon mass fraction first causes a drop in /cr at low 
temperatures as C/O approaches 1 (X('^C) x 2.0, dotted line), because more oxygen atoms get bound in CO and less other molecular 
opacity carriers can be formed. When C/O rises beyond 1 this trend is reversed and the mean opacity increases again due to the 
formation of carbon-bearing molecules. At higher temperatures the opacity is growing monotonically due to CO, CN and atomic 
carbon. See also Fig. [8] All curves have been smoothed using cubic splines. 



that are important to the Rosseland opacity are dissociated at 
these temperatures. The partial pressures of these molecules rise 
monotonically with the carbon abundance and thus the opac- 
ity at high log T increases in the same manner At intermediate 
temperature (around \ogT = 3.4) where many different opac- 
ity sources contribute, the aforementioned mechanisms compete 
and the behaviour of /cr is a more complex function of C/O. 

Briefly summarised, changes in the chemistry alter the 
Rosseland mean opacity, where variations in the C/O ratio 
has the most pronounced effect. An oxygen-rich composition 
(C/O < 1) results in a different group of molecules account- 
ing for the opacity than in the carbon-rich case (C/O > 1). 
Furthermore, carbon-bearing molecules have different spectral 
appearances from the oxygen-bearing ones and thus cause kr to 
have a different functional behaviour. iHelling et al.l (i2000) pro- 
vided examples of monochromatic absorption coefficients for 
carbon-bearing molecules. The only two molecules that con- 
tribute significantly for either chemistry are CO and CN. 

As emphasised earlier, the presented data are primarily rele- 
vant to the envelopes of evolved low mass stars. Chemical com- 
position variations due to the TDU acting in AGB stars concern 
mostly the enrichment in carbon. However, the dredged-up car- 
bon in the envelope can be inputted back into the CN cycle, 
which partly converts '^C to ''^N (see Introduction). By varying 
the abundance of nitrogen (more precisely '"^N), we add a fur- 
ther dimension to our data tables. These alterations have more 
direct consequences on the behaviour of the Rosseland mean in 
the sense that an increase in nitrogen always causes an increase 
in tiie opacity (in contrast to an increase in carbon, which can 
lower A-R for a certain parameter range; see above). Nitrogen is 
present only within two molecules being considered here, i. e. 
CN and HCN, and can directly influence the Rosseland opac- 
ity only via these compounds. Other molecules containing nitro- 
gen indirectly affect the opacity by altering the molecular par- 
tial pressures in chemical equilibrium. In the oxygen-rich case, 
the effect of an increase in ''*N is relatively moderate, because 
the partial pressure of CN is in general low. In the carbon-rich 
case, the abundance of CN is however much higher, HCN is also 



present in significant amounts, and, as a consequence, the opac- 
ity can increase considerably. These properties of kr are illus- 
trated in Fig. [3]for the two enhancement factors of ''^N. On top, 
we present results for two different metallicities without any car- 
bon enhancement, where the increase in opacity is due to CN 
only. In the respective bottom panels, the carbon abundance has 
been enhanced to its maximum value and HCN causes a con- 
siderable rise in kr at lower temperatures. In the high temper- 
ature range, some minor contributions from atomic nitrogen to 
the opacity are evident. 

The results discussed above are condensed into the form of 
14 separate files, one for each metallicity0Each file consists of a 
header indicating the abundances used, the initial metallicity, the 
initial mass fractions for '^C, '"^N, and the alpha elements, and a 
look-up table for the true data block. The final file consists of 63 
rectangular data arrays, where log kr is tabulated as a function 
of logT and log/?. The tables are ordered such that the mass 
fraction X('2C) varies the most rapidly followed by the hydrogen 
mass fraction and X(''*N). For future compatibility, a data field 
for the alpha element enhancement factor was introduced into 
the look-up table. 

5.1. Comparison with otiierdata 

We compare our tables based on a scaled solar metal mixture 
with data from F05 based on the same abundances as in this 
work. A direct comparison with AF94 is not possible b ecause 
there are, of course, no tables based on the lLoddersI (12003b abun- 
dances. We refer to F05 for a comparison of AF94 and F05. In 
the figures, we always depict data from our database as a-coma, 
while the respective comparison values are labelled kr. Despite 
the numerous differences between the COMA and F05 approach, 
we find reasonable agreement between both sets of data. For the 



^ The tables are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (13Q. 79. 128. 5) 
or via http : / /cdsweb . u-strasbg . f r/cg i-bin/qcat? J /A+a7] 
Alternatively, the files will be provided on request. 
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Fig. 3. Effects of an increased nitrogen abundance on the Rosseland mean opacity at various log/? values relative to the case 
without any nitrogen enhancement. The filled and empty symbols refer to intermediate and maximum enhancement factors of ''*N, 
respectively. For Z = 0.02 (left panels), these are 1 .5 and 2, whereas for Z = 0.00001 (right panels), we have 40 and 1500. On the 
top, we show the cases for the oxygen-rich case (default carbon abundance). Here the increase in opacity is due to the CN molecule. 
The bottom panels refer to maximum enhanced '^C (i. e. the carbon-rich case), where HCN also contributes to the mean opacity at 
lower temperatures. Note the different scales on the y-axis. 



case shown in Fig. |4] (Z = 0.02), the difference between the 
COMA and F05 values does not exceed 15 per cent for tem- 
peratures as low as log T - 3.5. The discrepancies at lower tem- 
peratures are higher (up to 35 per cent) and can in fact be as- 
cribed to several things. First and foremost, the use of different 
sets of molecular data in the calculations (cf. our Table |2] and 
their Tables 3 and 4) produces a deviation in the resulting mean 
opacity coefficients. Second, we adopt a microturbulent velocity 
of 2.5 km s"', while F05 use 2.0 km s"'. The choices for this pa- 
rameter are (within a certain range that is found for atmospheres 
of low mass giants) somewhat arbitrary and cause perceptible 
changes in /cr, especially at lower temperatures. Third, F05 use 
a denser wavelength grid for the evaluation of a-r. We discuss 
these issues in more detail in Sect. 15.31 From a comparison of 
Fig.|4]with Figs. [TOl (showing a comparable order of magnitude 
of the deviations) and[TT] it is, however, clear that the numerous 
differences in the physical input data are responsible for the ma- 
jor part of the discrepancies. The resolution and microturbulent 
velocity influence /cr not quite as much. The large deviations 
in the data at the lowest temperatures are due to grain opacity 
that we do not take into account in our calculations, but dust is 
usually not formed under equilibrium conditions (as assumed by 



F05, see Introduction). Moreover, F05 adopted a finer grid in 
log T below 3.5. For the oxygen-rich case, a cubic spline inter- 
polation (see Fig.m dotted lines) on the coarser grid we adopted 
(and also used by AF94) provides reasonably accurate values. 

The comparison with high temperature data such as that 
from OPAL or OP is limited to the temperature regions where 
the tables overlap. Moreover, it is this region where a transi- 
tion between low and high temperature opacities has to be made 
for applications covering a wide temperature range. OP data 
stretch down to log T - 3.5, whereas the OPAL tables end at 
log T - 3.75. The comparison for a standard scaled solar com- 
position in Fig.|5]shows a growing deviation for lower tempera- 
tures because both OPAL and OP do not include molecular ab- 
sorbers (except H2). This plot indicates that in the region be- 
tween log r = 3.8 and the high temperature end of the COMA 
data, a smooth transition to high temperature data is possible. 
Again, from the dimension of the differences, we conclude that 
these are due to different physical input data rather than other 
parameters (see Sect. 15.31 ). To assess which temperature region 
renders itself to such a crossover we plot the logarithmic dif- 
ference between our opacities and those of OPAL and OP, re- 
spectively, as a function of \ogR in Fig. |6] The deviations are 
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Fig. 8. Evolution of the Rosseland mean as a function of the C/O ratio for different metallicities (Z = 0.02 and Z = 0.00001) and a 
representative value of log/? = -1.5. To obtain opacities at enhanced carbon abundances in-between grid points we recommend a 
linear interpolation scheme in log /<-r and logX('^C) as indicated by the lines. The largest relative errors have to be expected at low 
metallicities and temperatures between X('^C) x2.2 and the successive enhancement factor. For a description of the mechanism that 
makes /cr drop to a minimum value at C/O = 1 and then rise again we refer to Sect. |5j 



in general moderate and the closest agreement is found around 
temperatures of 10,000 K. A few comments about the differences 
are required here. For a metallicity of Z = 0.02, the values from 
COMA are systematically higher than both OPAL and OP. This 
discrepancy almost vanishes at the lowest metallicity considered 
(Z = 0.00001, Fig.|7J. It must therefore be related to the metal- 
licity, which leaves either the metal lines or the differing contin- 
uum opacity as a reason. As for the first possibility, OPAL and 
OP use a restricted set of the most abundant elements, while we 
use all data that is contained in VALD. Whether this can cause 
the differences is a question to be answered by further analysis. 
Moreover, either OPAL and OP data do not include line broaden- 
ing due to microturbulence, although this alone cannot account 
for the difference in a-r. 

Beside the constant offset, we recognise a growing discrep- 
ancy at increasing log R, which is supposedly to be linked to the 
pressure broadening of the atomic lines. As mentioned earlier, 
the damping constants are those from VALD. Additionally, the 
adopted tabulated hydrogen line profiles could play a role since 
hydrogen lines contribute significantly to the Rosseland mean at 
high temperatures (cf. e. g. AF94). A deeper investigation of this 
issue is still to be completed, although the differences remain on 
an overall satisfactorily small scale. 

For a mixture enriched in carbon and nitrogen, we also 
performed tests. We switched to the lowest metallicity in our 
database where the effects should be most pronounced due to the 
high enhancement factors. Here we restrict ourselves to a com- 
parison with OP data since these can be produced quite easily 
using the software contained in the OPCD (version 3.3). Based 
on the same mass fractions that we used in COMA, we gen- 
erated opacity tables with the aforementioned program set and 
compared them with ours. The results are shown in Fig. [T] We 
refer to the discussion in the above paragraph concerning the dif- 
ferences but we emphasise again that the quantity logA-R/zc-coMA 
measuring the differences in the logarithmic opacity in dex re- 
mains within a reasonable range, i. e. A < 0.05 dex. In these 
plots, it is evident that the transition from the COMA coefficients 
to OP data should occur around a temperature of log T - 4.0. 
From the work of F05 who provided more details about the 
relation between the low and high temperature opacities, and 



ISeaton & Badnelll (12004 ?) who completed an in-depth compari- 
son between OP and OPAL, we conclude that the above state- 
ment also holds for OPAL data. 

5.2. Interpolation 

The interface between the opacity tables and the various applica- 
tion codes in which they are used are the interpolation routines. 
Concerning the temperature and log R, it has become sort of a 
standard to interpolate in these dimensions using cubic splines 
(e. g. ISeatonlli993[ 1 19961) or quadratic fits (in the Fortran sub- 
routines from Arnold I. BoothroycQ)- The results are quite satis- 
factory, but we want to emphasise that, in the carbon-rich case, 
problems can occur At low values of log R and a high carbon 
enhancement, a cubic spline interpolation in the log T dimen- 
sion might overshoot and produce spurious results. We strongly 
advise always checking separately the quality of the fit for each 
table (or relevant parts thereof) used. 

The problem now is how to account for the element enhance- 
ments. As outlined in Sect.|5] the special role of the C/O ratio the 
parameter range at C/O = 1 at low temperatures, and kr is not a 
continuous function of the carbon content at this point (Fig.[8]l. 
To resolve the sharp turnaround in the Rosseland mean, we re- 
quire some grid points close to C/0 = 1 . Overall, the number of 
enhancement factors is too low and the grid too coarse to apply 
any other interpolation scheme than a linear one. As shown in 
Fig. [8] for the solar case and at low metallicity, linear interpo- 
lation in log/fR and logX('-C) delivers quite gratifying results 
beyond a certain temperature, where molecules cease to play an 
important role in determining the value of the Rosseland opacity. 
The lower the temperature becomes, the sharper the turnaround 
in the functional behaviour of /cr. Due to the sudden drop in 
opacity when the amount of carbon and oxygen atoms become 
approximately equal, linear interpolation misses out a certain 
fraction of information. At high metallicities the situation is not 
so serious, although the case shown in the right panel of Fig. |8] 
(Z = 0.00001) reveals this shortcoming clearly. Once the ele- 
ment mixture is carbon-rich (log C/O > 0), the opacity first in- 
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Fig. 4. Comparison of COMA (filled circles) and F05 (empty tri- 
angles) values at a metallicity of Z = 0.02. The top panel shows 
absolute values, whereas in the bottom panel the differences be- 
tween F05 and COMA values are indicated on a relative scale. 
Here and in the following figures, /ccoma refers to data contained 
in our database, while the respective comparison values are al- 
ways labelled k^. The overall agreement down to logT = 3.5 
is gratifying. The growing discrepancies towards lower temper- 
atures are most probably due to the deviating set of molecular 
data in the respective calculations. The steep increase in in 
the F05 data below log T = 3.25 is due to dust grains that are not 
accounted for in our data. Below log T - 3.5, F05 use a smaller 
temperature spacing than we do. The regions between the grid 
points can be reasonably well reconstructed by a cubic spline 
interpolation in log T (dotted lines in top panel). In the bottom 
panel, we compare values at the COMA grid points only. 



creases sharply but flattens at high carbon enhancement values. 
Due to the construction of our tables (see Sect. |4|, the spacing 
of the enhancement factors in the carbon-rich regime increases 
at lower metallicities. Between the X{^^C) x 2.2 and the suc- 
cessive enhancement factor, an additional grid point would be 
favourable. For the case of nitrogen enhancement, linear inter- 
polation in both log /fR and logX('"*N) is a good approximation, 
because the relation between the nitrogen content and the opac- 
ity has a simpler behaviour. 

5.3. Sources of uncertainties 

The definition of the Rosseland mean opacity in Eq. [T] leaves 
only some ambiguity about how to evaluate this quantity in terms 
of numerical methods. However, the considerable uncertainties 
in published opacity coeflicients originate in data entering the 
calculations. In the case of low temperature opacities, there is, 
in particular, a good amount of physical data of different quality 
that must be combined into one quantity. The summary in the 
following paragraphs is not exhaustive but discusses the accu- 
racy of the data presented here and elsewhere. 
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Fig. 5. Comparison of COMA values with OP and OPAL for 
Z = 0.02. The top panel shows absolute values, whereas in the 
bottom panel the differences in the overlapping region are again 
indicated on a relative scale. The high temperature data do not 
include the relevant molecular opacity sources. Therefore, the 
discrepancies increase the lower the temperature becomes. OP 
data range down to log T - 3.5, whereas the OPAL tables end at 
log T - 3.75. The transition between low and high temperature 
data can be done close to the upper temperature border of the 
COMA data, where agreement between the values shown here is 
fairly good. See also Figs.|6]and|2] 
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Fig. 6. Comparison of COMA values with OP and OPAL. Our 
opacity coefficients are systematically higher for a metallicity of 
Z = 0.02. This constant offset vanishes for lower metallicities 
(see Fig.|7]i, which implies that atomic opacities (lines and con- 
tinuum) of the metals are the cause. For high values of log/? 
where the pressure broadening of the spectral lines becomes 
important, we find increasingly divergent results. The plots ap- 
pear to imply that the low and high temperature data should be 
merged around log T - 4.0 (solid line). 
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Fig. 7. Comparison of COMA with OP at a low metallicity 
(Z - 0.00001) with scaled solar abundances of carbon and nitro- 
gen, and with C and N enhanced to the maximum values (which 
results in a higher metallicity). The data agree quite well, the 
discrepancies remain within 5 per cent around 10,000 K (solid 
line). Again, at high log/? the deviations are higher, obviously 
due to a deviant treatment of the pressure broadening of spectral 
lines. From a comparison between OP and OPAL (see text) we 
draw the same conclusions for OPAL data. 



5.3.1 . Adopting different solar element abundances 

Hitherto, it has been emphasised that low temperature Rosseland 
opacities are to a large extent determined firstly by the total 
metallicity of the element mixture Z and secondly by the C/O 
ratio. The individual element abundances play a role as well 
but do not change the opacity on an order-of-magnitude scale 
when the aforementioned parameters are kept fixed. Concerning 
the oxyge n-rich case, we ref er to a discussion of this topic 
given by iDotter et"al] ( |2007|) . For the carbon-rich case, we 
exemplarily calculated opacity coe fficients starting from so- 
lar element abundances other than iLodde rs' 1*2003 1 L03), that 
is iGrevesse & Sauvall (Il998l GS98) and [Orevesse etal] (l2007l 
GAS07). The results of our comparison are shown in Fig. |9] 
(for Z = 0.02 at logR = -3.0). The solar abundances given by 
GAS07 are very similar to those of L03, and it is therefore unsur- 
prising to uncover virtually identical opacity coefficients for our 
test case. The situation is different when we consider the GS98 
abundances, which provide higher values for C, N, and O than 
L03 and GAS07. These elements also make up a higher fraction 
of Z than in the other cases. In turn, when the metals are scaled 
to obtain Z = 0.02, metals apart from C, N, and O, are present 
in lower amounts than in the L03 case, and thus contribute a 
smaller fraction to the opacity at high and intermediate temper- 
atures (beyond log T = 3.5). At the lowest temperatures, more 
carbon-bearing molecules, such as C3 and C2H2, are likely to 
form and produce a higher value of /cr, partially compensating 
for the lower atomic opacity contribution at intermediate tem- 
peratures. The sensitivity of our results to the adopted starting 
abundances is limited. The size of the differences with respect to 
the standard case is similar to other uncertainties discussed here. 
Hence, our data can be used to approximate the Rosseland opac- 
ity coefficients for a different set of scaled solar abundances, as 
long as the CNO abundances do not deviate significantly from 
the values given in L03. 
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Fig. 9. Changes in the Rosseland opacity when using 
different sets of s olar element abu ndances, for i n stance 
IGrevesse & Sauvall ( Il998i GS98) and iGrevesse et all ( l2007l 
GAS07). Data shown here are for the carbon-rich case and 
Z - 0.02 at log/? = -3.0. The C/O rati os were s et to match 
the values contained in our database using iLoddersl (i2003. L03) 
abundances resulting from ^('^C) x 2.2,3.0,5.0. GAS07 is in 
many respects very similar to L03 (e. g. regarding the C, N and 
O abundances and the share of these elements in Z) and thus 
results in almost identical opacity coefficients. The abundances 
given by GS98 deviate much more from L03. The differences 
become manifest in the opacities at high temperatures, and also 
at low log T by means of the chemical equilibrium. Other uncer- 
tainties discussed in the text have a comparable order of magni- 
tude. 
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Fig. 10. Uncertainties in the molecular data that affect the mean 
opacity. Top panel: Using a different line list for water which is 
the main source of opacity at low temperatures in the oxygen- 
rich case causes considerable changes in /cr. As a test we substi- 
tuted the BT2 water line list by the one from SCAN and plot the 
relative changes compared to the default setup. Bottom panel: 
For the database described in this paper we scale the C2 line 
strengths in a certain wavelength region (see text). Not apply- 
ing this correction results in notable differences in the Rosseland 
mean opacity of carbon-rich mixtures. (X = 0.7, Z = 0.02, no 
and maximum carbon enhancement in the top and bottom panel, 
respectively.) 
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5.3.2. Uncertainties in molecular data 

We first consider the molecular line data. For many molecules, 
there is more than one line list available. The Rosseland mean as 
a global quantity is insensitive to the precision of line positions 
so long as the overall opacity distribution is reproduced well. 
However, there are cases where data from different sources result 
in altered overall opacity coefficients. As an example, we con- 
sider the contribution of water, the major low temperature opac- 
ity source in the oxygen-rich case, to kr. In this work we us e the 
recently published BT2 water line list jBarber et alj|2006l) . An 
alternative would have been the SCAN database line list from 
[j0rgensen et al. (2001). For the solar metallicity case, we cal- 
culated a table utilising this list and illustrate the results in the 
top panel of Fig.fTO] The discrepancy in the resulting values are 
as high as 30 per cent. Pronounced differences between the two 
lists lie in the region where the weighting function in the defini- 
tion of the Rosseland mean has its maximum, and we ascribe the 
deviating values of kr to this fact. 

For the uncertain line data in the carbon-rich case, we men- 
tion th e modifications to the C2 line data from Ouerci et alj 
(II974I) . To reproduce carbon star spectra. lLoidl et al.l (1200 lb pro- 
posed a scaling of the g f values in the infrared region (suggested 
by fj0rgensen 1997) based on a comparison with other line lists. 
More precisely, they scaled the line strengths by a factor of 
0.1 beyond 1.5 yum and left them unchanged below 1.15 yum. In- 
between, they assumed a linear transition. We adopt this method 
for the calculation of our opacity tables. By not applying this 
modification to the line strengths, we would have caused an in- 
crease of kr of roughly 25 per cent at Z = 0.02 with maximum 
enhanced carbon (Fig.[TOl bottom panel). The error in these data 
will have a more significant effect at low metallicities, where one 
expects a higher enrichment in carbon. From the calculation of 
mean opacities, we observe a clear need for new and improved 
C2 fine data. 

Beside the problems with existing data there are also 
molecules so far unconsidered that are suspected of providing 
non-negligible contributions to the opacity. The prime example 
is C2H, which could be an important opacity source in carbon 
stars, altho ugh to date no Une data has existed for this molecule 
(we refer to iGustaf ssonll 1 9951 for an overview). 

Another decisive set of input parameters are the chemi- 
cal equilibrium constants usually depicted by Kp. Each con- 
stant is in fact a temperature-dependent function setting the par- 
tial pressure of a molecule in relation to the product of the 
partia l pressures of t he mo lecule's constituents (cf. e. g. iTsuiil 
Il973h . iHelUng et all ( |2000|) pointed out that the literature val- 
ues for equilibrium constants from different sources could differ 
strongly at low temperatures. The critical point here is that one 
has not only to pay attention to the main opacity carriers but also 
to less abunda nt molecules compet ing with them for the same 
atomic species. iHelUng et alj ( 1200 0) referred to TiO and Ti02 as 
a examples but also reported other molecules for which order-of- 
magnitude differences in the partial pressures were found using 
different sets of Kp data. The data that we use is documented in 
Sect. [3] 

The above examples underline that accurate molecular line 
data is not only desirable for high resolution applications but also 
of importance to the calculation of mean opacities. In general, all 
data used in calculating the Rosseland mean, whether line data or 
other accompanying data like partition functions or equilibrium 
constants and continuum sources, must always undergo critical 
evaluation. 
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Fig. 11. Influence of variations in resolution or microturbulence 
on the mean opacity. Top panel: The calculation of kr requires 
to know the monochromatic opacity at a large number of wave- 
length points. A high spectral resolution assures convergence 
but induces high computational costs. We use a lower resolution 
than F05 but the changes in the opacity coefficients are rather 
moderate especially when compared to other sources of uncer- 
tainty. Bottom panel: We adopted a microturbulent velocity of 
^ = 2.5 km s ' throughout this work. Changing this value to 
^ = 2.0 km s"' like in F05 causes kr to drop. (X = 0.7, Z = 0.02, 
no carbon enhancement.) 



5.3.3. Numerics and parameters 

Apart from the imprecision due to the physical input data, there 
are other factors influencing the emerging opacity coefficients. 
For instance, an error source exists in the wavelength grid for 
which the opacities are calculated and the integration completed 
in deriving the opacity mean. Compared to F05, we use a con- 
siderably lower spectral resolution. To assess the uncertainties 
due to this difference, we simulated the resolution of F05, recal- 
culated one of our tables, and compared our results to those for 
the original case. The differences found were relatively small as 
shown in Fig. [TTI (upper panel). Since the error is low compared 
to other effects described above, we propose that use of a lower 
resolution is justifiable because it would reduce considerably the 
amount of CPU time. 

On the other hand, additional physical parameters enter the 
calculation of a-r, such as the microturbulent velocity ^, which 
influences the width of the line profiles. The spectral lines 
are broadened according to the adopted value for ^, which is 
somewhat arbitrary. Throughout this work, we used a value of 
^ = 2.5 km s ' for the generation of our data. Results from 
previous works on spectra of late-type sta rs (e. g. lAringer et al.l 



2002; Gautschy-Loidl et al. 2004; Lebzelt er et all l2008l) have 



shown that this is a reasonable assumption. In the work of F05, 
however, ^ was set to equal 2.0kms '. Both options are well 
within the range of valu es found for AGB star atmospheres (e. g. 
ISmith & Lamber3[T990l) . In Fig. [m (lower panel), we show the 
results of a test using the F05 value. Since the spectral lines pos- 
sess a smaller equivalent width at a reduced value of ^, the mean 
opacity is lower than for the COMA default case. 
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5.3.4. Application of tine data 

Beyond generating tables, there is further possibilities for in- 
troducing uncertainties while applying the data. First, there 
is the technical problem of interpolating the tabulated values. 
Compared to previously available data, the situation is worse be- 
cause there are two more dimensions along which to interpolate, 
that is the varying amount of carbon and nitrogen. However, on 
the basis of the above discussion, it is unlikely that too sophisti- 
cated interpolation algorithms produce improved accuracy. This 
is, however, a problem that can in principle be solved by increas- 
ing the amount of computer power 

Far more worrying and the largest error source of all is po- 
tential misapplication of the data. Strictly speaking, the scope of 
the tables containing Rosseland mean opacity coefficients are re- 
gions where the diffusion approximation for the radiative trans- 
fer is fulfilled. In terms of the optical depth, this means r » 1 
for all wavelengths. One of the main applications of our data 
will be the outermost parts of an AGB star evolution model. The 
outer boundary condition is usually set somewhere in the atmo- 
sphere (logT < 3.6), where by definition t < 1. In some situa- 
tions, the Rosseland mean might still be a good approximation 
for evaluating the radiative energy transport. However, in gen- 
eral it is necessary to use a non-grey radiative transfer method 
because, due to the molecular absorbers, the spectral energy dis- 
tri bution is strongly w avelength-dependent. We refer to the work 
of lHofner et al.l (l2003h . who demonstrated the shortcomings of a 
grey treatm ent of the radiative transfer f or dynamical model at- 
mospheres. [Harris & Lvnas-GrevI (l2007i) investigated the effect 
of non-grey surface boundary conditions on the evolution of low 
mass stars and reported noticeable changes to RGB evolution 
tracks. We thus want to emphasise that our mean opacity tables 
are meant to provide an interim solution until modelling of non- 
grey radiative transfer in stellar evolution calculations becomes 
feasible. 

6. Conclusions 

We have presented a grid of low temperature Rosseland mean 
opacity tables that take into account variations in the single 
element abundances of carbon and nitrogen. By gradually en- 
hancing the carbon content of a metal mixture, the molecular 
contribution to opacity changes significantly due to the altered 
chemistry. Akeady within a certain regime (i. e. oxygen-rich or 
carbon-rich), the relative amount of carbon to oxygen has pro- 
nounced effects on kr. More distinctive, however, is the com- 
parison between oxygen-rich and carbon-rich regimes. Different 
molecules serve as opacity sources in either case and thus result 
in a qualitatively and quantitatively different Rosseland mean 
opacity as a function of temperature and density. Changes in the 
nitrogen abundance also alter the opacity coefficients via certain 
nitrogen-bearing molecules. 

The tables are designed such that an incorporation into exist- 
ing codes that utilised AF94 or F05 data should be straightfor- 
ward. Our data cover a wide metallicity range, and the overabun- 
dances of carbon and nitrogen are adjusted in each case. We are 
confident that with these data we provide a tool to simulate the 
final phases in the evolution of low or intermediate mass stars 
in more detail. Our data include the effects of the ongoing nu- 
cleosynthesis and mixing events in AGB stars in terms of th e 
opacity. As shown in previous papers (e. g. lCristallo et al.ll2007h . 
the incorporation of our tables into stellar evolution codes alters 
the physical properties of the stellar models. Once the star be- 
comes carbon-rich, molecules form that are more opaque than 



those in an oxygen-rich regime. This in turn results in a steeper 
temperature gradient. A consequence is, for instance, a decrease 
in the effective temperature in stellar evolution models. The stel- 
lar radius increases, and the average mass- loss rate increa s es and 
erodes the envelope mass at a faster rate. Straniero et alj (|2003|) 
showed that a change in the envelope mass (as well as a change 
in the core mass) affect fundamental properties of AGB stars, 
e. g. the strength of the thermal pulses and the total amount of the 
mass dredged up. It will also be interesting to see how different 
mass-loss prescriptions interact with the newly calculated opac- 
ity coefficients, since these issues are physically closely coupled. 

In the future, we plan to extend our tables to contain data 
about the enrichment in alpha elements. We must emphasise, 
however, that the data provided in the course of the current and 
future work must be seen as a transitional solution to the treat- 
ment of molecular opacity in AGB star envelopes and atmo- 
spheres. Due to the band structure of molecular absorption, mean 
opacities will yield inaccurate results. The past results from 
static and dynamical model atmosphere calculations demon- 
strate the importance of a frequency-dependent radiative trans- 
fer. Our data promise to bridge the gap until these methods are 
employed in stellar evolution models. 

Finally, we emphasise that for forthcoming extensions of this 
database, it would be desirable to obtain extensive response from 
the community. Comments and criticisms that can lead to an im- 
provement in the quality of the data are highly welcome. 
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